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Dynamical quasi-localization of matter fields in extra dimension is an elegant mechanism to gen- 
erate hierarchical 4-dimensional Yukawa couplings. We point out that a bulk matter field whose 
zero mode is quasi-localized can give a large Kaluza-Klein threshold correction to low energy gauge 
couplings, which is generically of the order of ln(Yukawa)/87r 2 , so it can significantly affect gauge 
coupling unification. We compute such threshold corrections in generic 5-dimensional theories com- 
pactified on S 1 /Z2 x Z2, and apply the result to grand unified theories on orbifold generating small 
Yukawa couplings through quasi-localization. 

I. INTRODUCTION 



It has been noticed that extra dimension can provide an elegant mechanism to generate hierarchical Yukawa 
couplings [1-14]. The quark and lepton fields can be quasi-localized in extra dimension in a natural manner, and 
^sO • then their 4-dimensional (4D) Yukawa couplings are determined by the wavefunction overlap factor e - MlrR where 
M is a combination of mass parameters in higher dimensional theory and R is the length of extra dimension. This 
V/ allows that hierarchical Yukawa couplings are obtained from fundamental mass parameters having the same order of 
magnitude. Extra dimension has been known to be useful also for constructing a natural model of gauge unification 
[15-25]. Supersymmetric 4D grand unified theories (GUTs) successfully accomodate gauge coupling unification, 
£C} < however suffer from some difficulties such as the doublet-triplet splitting problem and the proton decay problem. 
These problems can be elegantly solved in GUT on orbifold with gauge symmetry broken by boundary conditions 
[15, 18] . It is straightforward to implement the idea of dynamical quasi-localization in orbifold GUT to get hierarchical 
£^ I Yukawa couplings as well as successful gauge unification [7, 9]. 

In any GUT, heavy particle threshold effects at GUT symmetry breaking scale should be taken into account for a 
precision analysis of low energy gauge couplings g 2 a . In conventional 4D GUT, those threshold corrections to l/g% are 
generically of the order of 1/87T 2 and thus not so important, unless the model contains a large number of superheavy 
particles which become massive as a consequence of GUT symmetry breaking or some of superheavy masses are 
hierarchically different from each other [26]. As was pointed out a long time ago, higher dimensional field theory 
and/or string theory contain (infinitely) many Kaluza-Klein (KK) or stringy modes, so can have a sizable threshold 
correction [27]. Therefore it is essential to include stringy and/or KK threshold correction in the precision analysis 
of low energy couplings in string and/or higher dimensional field theories. In this paper, we wish to examine the KK 
threshold corrections in generic 5D orbifold field theories in which hierarchical 4D Yukawa couplings are generated 
through quasi-localization. As we will see, the KK threshold corrections to \jg\ in such models are generically of the 
order of ln(Yukawa)/87r 2 , so can significantly affect gauge coupling unification. The outline of this paper is as follows. 
In Sec. II, we discuss the KK spectrum and zero mode wave functions of generic scalar and spinor fields in 5D theory 
compactified on S 1 /Z2 x Z' 2 . In Sec. Ill, we discuss the 4D Yukawa couplings of those quasi-localized scalar and 
fermion zero modes. In Sec. IV, we compute KK threshold corrections in generic 5D theories on S 1 /Z2 x ^. In Sec. 
V, we use these results to derive the KK correction to the predicted value of the low energy QCD coupling constant 
in 5D GUT on orbifold. We then construct a class of supersymmetric 5D orbifold GUTs which generate hierarchical 
Yukawa couplings through dynamical quasi-localization while keeping successful gauge coupling unification. Sec. VI 
is the conclusion. 
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II. KALUZA-KLEIN ANALYSIS AND QUASI-LOCALIZED ZERO MODES 

In this section, we analyze the KK wave functions and spectrums of scalar and spinor fields in 5D theory compactificd 
on S 1 /Z 2 x Z' 2 . Our major concern is the dynamical quasi-localization of zero mode wavefunctions which would result 
in hierarchical 4D Yukawa couplings. We also consider the full KK spectrums which will be relevant for the discussion 
of KK threshold corrections to low energy gauge couplings. Throughout this paper, we will use the convention for 
S 1 jZ 2 x Z' 2 which is given by the following transformation of the fifth spacetime coordinate y = y + AnR: 

Z 2 :y^-y, Z 2 : y + ttR -> -y + ttR . (1) 

so the fundamental domain of S 1 /Z 2 x Z 2 corresponds to < y < nR. 
Let us first consider a 5D complex scalar field with an action: 



/ 



D M ^'D M ^-'* + ( m 2 zz , + 2^ zz ,5(y) - 2^ zz ,5(y - nR) ) <t> zz ' 4> zz '* j , (2) 



and the orbifold boundary conditions: 

4>"' \-y) = z<t> zz ' \y) , 

<j> zz \-y + nR)=z'$ zz '{y + nR), 

where z — ±1 and z' — ±1. The equation of motion for each KK mode of this scalar field is given by 



(dl + K z ') 2 ) ft*' = 2 [» zz ,8(y) - n' zz ,6{y nR) ] , (3) 

where 

(oof) 2 = (M„"') 2 -m\ z , (4) 

for M zz denoting the 4D mass of <p^ ■ Upon ignoring unimportant normalization factor, one easily finds the following 
KK wave functions on the fundamental domain: 

4>++ = cos {uj+ + y) + sin (LU++y) , 

4>+~ =sin(w+-(j/-7rii)) , 
<j>- + = sin {u)~ + y) , 

</>-- = sm(w--y) , (5) 
where cj^ z 0) are either real or pure imaginary constants determined by 



tan (uj+ + ttR) 


K + ) 




tan (cj+ _ 7ri?) 


M+- 




tan (w~ + 7ri?) 






tan (u)~~irR) 


= 0. 





(6) 

Once oj^ z ' are determined, the KK mass spectrum of <j> zz ' can be determined by the relation 

(M-') 2 = mL, + (<') 2 . (7) 



Note that a pure imaginary uj^ z gives (m zz ) smaller than m 2 zz , , while a real u>^ z gives (m zz ) bigger than m 2 zz , . 
Note also that there can be only one or no imaginary lo zz from (6). 
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Since we are interested in the shape of zero mode wavefunction, let us consider the conditions for <f> zz to have a 
zero mode, i.e Mq z = 0, independently of the value of the orbifold radius R. Obviously only </> ++ can have such a 
zero mode only when its bulk and brane masses are tuned to satisfy 

M++ = »++ > m ++ = /4+ , (8) 

which are indeed satisfied in supersymmetric theories [5, 28]. The resulting zero mode (on the fundamental domain) 
is given by 

</>++ = exp(/i ++ y) , (9) 

so it is quasi-localized at y = if ^ ++ < 0, while at y = irR if > 0. In a more general case with = 
but m 2 lr+ 7^ the KK spectrums of <p ++ are given by 

(M++) 2 = m 2 ++ -4 +) 

(M+ + ) 2 =ml + + ^, (10) 

for n being a positive integer. 

For </> + ~, if ttR < — 1, there can be a mode lighter than : 

(M+-) 2 = m 2 + _-e 2 , (11) 

where £ is a real constant determined by 

tanh(£7ri?) = . 

|M+- 1 

In the limit 7r_R <C — 1, Mq is determined to be 

(M+-) 2 = (ro*_ +4/4_ ( e - 2 l"+-l-« + 0(e- 4 l' t +-l 7rfl )) . (12) 

So in supersymmetric case in which m+_ = the mass of <j>Q _ is exponentially suppressed in the limit fi^ ttR <C 

— 1. On the other hand, all other KK modes of H are heavier than m + _. Similarly, if /j, ^irR > 1, h can have 

a lighter mode: 

<Po + = sinh(£'y) , 

(M -+y = m *_ + - e , 



where 



In the limit fj, ^ttR >• 1, we have 



tanh(£V.R) 



M-+ 



(M-+) 2 = (m 2 _ + - M 2 _ + ) 2 + 4 M 2 + ( e - 2 "-+** + 0( e -^-+^)) , 

yielding an exponentially small mass of 4>q + in supersymmetric case with /j, ^ttR 3> 1. 

Let us now consider a 5D spinor field with an action 

- J d 5 x [w zz ' ( 7 M 9 M + M zz <e(y)) V zz ' ] , (13) 

and the orbifold boundary condition 

* zz '(-y) = zj 5 * zz '(y), 

* zz '(-y + wR) = z'j 5 ^ zz '{y + nR) 7 
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where e(y) = —e(—y) = e(y + 2ttR) = 1 on the covering space of S 1 /Z2 x Z' 2 . Note that in order for the action to be 
invariant under 2 2 x Zj, 5D spinor can have only a kink type mass M zz >e(y). The equation of motion for each KK 
mode of ^ zz is given by 

(d 2 y + (ufrfni = -2M ZZ , [5{y) - S(y nR) } K*l , 

(dl + K/') 2 ) 2 = +1M ZZI [6(y) S(y nR)]K z R , (14) 
(^) 2 = (MD 2 -(M ZZ ') 2 (15) 



where 7 5 *l,a = ±*l,a and 



for M zz ' denoting the 4D mass of ty^' . When compared to the scalar equation of motion (3), this suggests that in 
supersymmetric case the bulk and brane scalar masses should satisfy 

M*z'=/4*'. ML' = m lz' =M\ Z ,. 
Obviously, has a left-handed zero mode for any value of A4++: 

*++(y)=exp(-M ++ y) (16) 

which is quasi-localized at either y = (if M++ > 0) or y = nR (if M++ < 0). All other left-handed modes of 
are paired-up with right-handed modes to get a 4D Dirac mass bigger than M++: 



n 2 



M++ = ^M\ + + '- 2 . (17) 

Similarly, has a right-handed zero mode for any value of M and also the massive modes with (M~~) 2 = 

M 2 __ + n 2 /R 2 . The wave function of zero mode is given by 

*ofl(y) = exp(X— y), (18) 

so the zero mode is quasi- localized at either y = (if A4 < 0) or y = irR (if > 0). 

For \D ,+_ , there is no zero mode. However if M+-nR > 1, there can be two modes lighter than , while all 

other modes are heavier than The wavefunctions of light modes are given by 

vf+ L -(y) = smh(k(y-nR)), (19) 
*0fl(y) = smh(fcy), (20) 

(21) 



where k is a real constant determined by 



tanh(fc7ri?) = — ^— . (22) 
M+- 



Obviously, one of these two light modes is localized at y = 0, while the other is localized at y = irR. In 4D viewpoint, 
these two modes are paired up to get a Dirac mass 

(M+-) 2 =M\_ -k 2 . (23) 

Note that in the limit M + -nR > 1, 

M+- w 2M+- cM-M+-ttR), (24) 

so can be arbitrarily light. Similarly, if A4 ^irR < —1, h can have a light Dirac mode with Mq + w 

2\M-+\cxp(M-+ttR). 
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The dynamical quasi-localization of light modes and also the shape of full KK spectrums offer a possibility that 
massless brane fields can be considered as a large mass limit of bulk fields 1 . To see this, let us take the limit 
\M. zz i | — > oo 2 . In this limit, the left-handed zero mode of \I>++ becomes a chiral brane fermion confined at y = 
(if M++ — > oo) or at y — ttR (if M++ — > — oo), while the right-handed zero mode of \f ( becomes a chiral brane 
fermion at y = (if A4— — > — oo) or at y = 7ri? (if .M__ — > oo). For the light modes and ^q^" become a 

massless brane fermion at y = and y = 7r_R, respectively, in the limit .M-i — ► oo. Similarly, ^ol" anc ^ ^oi?" become 

a massless brane fermion at j/ = 7ri? and y = 0, respectively, in the limit M. |_ — -> — oo. All other KK modes of , I ,ZZ 

are heavier than |.M ZZ '|, and thus are decoupled in the limit \A4 ZZ >\ — > oo. So any massless brane fermion can be 
considered as a bulk fermion in the large kink mass limit. 



III. YUKAWA HIERARCHY FROM QUASI-LOCALIZATION 

To discuss flavor hierarchy arising from quasi-localization, let us consider a generic 5D theory containing arbitrary 
number of scalar and fermion fields: 

S = -Jd 5 x [ Dm4>iD m 4>*i + (mjSu + 2 f i IJ 6(y) - 2 f i' IJ 6(y - nR)) fa fa, 

+ i^ A { 1 M D M +M A e{y))^A+C Y ] . (25) 

Here Ly stands for the Yukawa couplings between fa and * a and the orbifold boundary conditions are given by 

fa(-y) = zifay), fa(-y + nR) = z'jfa{y + nR) , 
^a(-v) = za!^a{v), ^ A {-y + TiR) = z' A ^ A {y + TTR). (26) 

Though the brane masses of scalar fields can have off-diagonal components in general, we will assume for simplicity 
that they are diagonal: /i/j = /i/<5/j and /i' 7J = [i'jSij. To get hierarchical 4D Yukawa couplings, we also assume 
that Yukawa couplings in 5D theory exist only at the fixed points, which is assured by 5D SUSY in supersymmetric 
theories. Then the most general form of Yukawa couplings can be written as 

A A' 

Cy = 6 (f) Vx ^ P ^ Q +S(y- nR) -^y^ ip x V> P i) Q +c.c, (27) 

where A denotes the cutoff scale of our 5D orbifold field theory, ip x — {tf>j or fa}, and ip p = {^(1 + 75)^ a or ^(1 + 
75) V S ,C A }. Note that \i pq and X' ipq are dimcnsionless parameters in our convention. 

Let 4n denote 5D scalar fields having a zero mode fai = exp(mjy), i.e. scalar fields with z% = z\ = 1 and 
Hi = [i\ = nii, and ^ p denote 5D spinor fields having a chiral zero mode tp np = cxp(—z p M p y), i.e. spinor fields with 
z p = z' p = ±1 . It is then straightforward to find that the 4D Yukawa couplings of canonocally normalized zero modes 
are given by 

Vr Vq = \J Z(mi)Z(-z p M p )Z(-z q M q ) X ipq 

+ y/z(-mi)Z(z p M P )Z(z q M q ) \' ipq (28) 

where 

*<»>-X?=anr (29) 

Obviously, y ipq can have very different values, depending upon the values of m, and M p , q , even when all of the 5D 
parameters X ipq and X' ipq have similar values. 

A simple way to get hierarchical Yukawa couplings through quasi-localization is to assume that all Yukawa couplings 
originate from a single fixed point, for instance from y = ttR. As a concrete example, let us consider the case that 



1 Here we are implicitly assuming that there is a 3-brane at each orbifold fixed point. 

2 This infinite kink mass should be considered as a large mass comparable to the cutoff scale of the orbifold field theory under consideration. 
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all Yukawa couplings originate from y = ttR and <pi is a brane field confined at y — ttR with Z(—m,i) w 1. Note that 
a brane scalar field at y = nR can be obtained from a bulk scalar field by taking the limit m = 0(A) while keeping 
fi = fi' = m. Then for z P:q M p . q < —1/R, we have 



Vi „ a ,,m^A Kn , (30) 



while for z p , q M p . q > 1/i?, 



y ipq „ ^4AV^ e -(^ P+ ^ 9 )^ (31) 

The physical interpretation of this result is simple. If z p , q M p . q > the corresponding zero modes are quasi- 

localized at y = 0, so the Yukawa couplings are exponentially suppressed as they originate from y = nR. On the 
other hand, for z p ^ q M. PA < — 1/R, the zero modes are localized at y = nR, so there is no exponential suppression in 
Yukawa couplings. 

If the fundamental theory at A is weakly coupled, simple dimensional analysis would suggest that the dimensionlcss 
Xi pq are of order unity or less. However a more interesting possibility is that the theory is strongly coupled at A 
[29, 30]. In 5D theories under consideration, the standard model gauge fields live in bulk spacetime, so their 4D 
couplings are given by \jg\ ~ nR/g§ a , where g| a denote the dimensionful 5D gauge couplings with mass dimension 
— 1. In order for our 5D theory to be a useful framework, it must be valid up to an energy scale significantly higher 
than the compactification scale 1/i?, i.e. A> If A were comparable to 1/R, we would not have to consider a 

5D theory as an intermediate step going from the fundamental theory at A to the 4D effective theory for low energy 
physics. We would rather go directly to the 4D effective theory from the fundamental theory. On the other hand, if 
A » 1/R as desired, we have g\ a w nRg 2 a ^> 1/A, implying that the fundamental theory at A is strongly coupled. As 
is well known, only in strongly coupled scenario, GUT on orbifold can provide a meaningful prediction for sin 2 9w or 
the QCD coupling constant at the weak scale. In such strongly-coupled scenario, dimensional analysis suggests 

K pq = 0{An). 

It is then straightforward to get hierarchical 4D Yukawa couplings ranging from the top quark Yukawa coupling y t w 1 
to the electron Yukawa coupling y e w 10~ 5 — 10~ 6 within the parameter range 

AnR = O(10 2 ) , \M p \ttR<7. 

Note that, in strongly coupled scenario, if the Higgs boson and the left and right-handed top quarks are all brane 
fields, the resulting top-quark Yukawa coupling would be too large, y t — 0(4n), so one needs to put some of those 
fields in bulk spacetime. 

In supersymmetric 5D theories, the fermion kink masses M p are related to the graviphoton gauging [32, 33]: 

Dm^p = d M ^ P + iM p e(y)B M «f p + ... , (32) 

where Bm denotes the graviphoton. This suggests that it is a plausible assumption that M p are quantized (in an 
appropriate unit) as the conventional gauge charges are quantized. If true, the resulting Yukawa couplings (31) would 
have the same form as those obtained from the Frogatt-Nielson mechanism which generates small Yukawa couplings 
using a spontaneously broken U(l) flavor symmetry [31]. In the next section, we will see that the KK threshold 
corrections to low energy gauge couplings are generically given by 

if the hierarchical 4D Yukawa couplings yi pq are generated by quasi-localization. Obviously then the KK threshold 
corrections can significantly affect the gauge coupling unification. 



IV. KALUZA-KLEIN THRESHOLD CORRECTION TO LOW ENERGY GAUGE COUPLING 

In this section, we discuss the 1-loop threshold correction to low energy gauge coupling in 5D orbifold field theory. 
The bare action of bulk gauge fields can be written as 

Sbarc = - J d 5 x ( Jjj- + ^5(y) + ^8(y nR)^J F aMN F a MN , (34) 
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and then the 4D gauge couplings at tree level are given by 



1 \ nR , 

— 2 I — — 2 I" K a + « a - 

9a J tree 55a 



(35) 



One-loop corrections to low energy couplings can be computed by summing the contributions from all KK modes. In 
generic 5D orbifold field theory, such computation yields [34, 35] 



skip) 



7a 

9la 247T3 
f \ 1 



+ 



1 

8^ 



A a (lnA,m,,u, fj,',M,R) + b a In ( — 

P 



9a 



8tt 2 



A a (]nA,m,fj,,fj,',M,R) + b a In 



(36) 



where j a are the coefficients of (UV-sensitive) linearly divergent corrections, A a stand for (UV-insensitive) logarithmi- 
cally divergent or finite threshold corrections due to massive KK modes, and b a are the standard 1-loop beta function 
coefficients due to zero modes. Here (l/g 2 )barc correspond to the uncalculable bare couplings of the model, while A a 
are unambiguously calculable within 5D orbifold field theory. Note that A a contain a piece linear in In A as well as 
a finite piece depending on the scalar and fermion mass parameters m, fi, // ,A4 and also on the orbifold radius R. 
Here we assume that the cutoff scale A is large enough compared to other mass parameters of the theory, thus ignore 
the part suppressed by an inverse power of A. In the above expression, the renormalization point p is assumed to be 
below the mass of the lightest massive KK mode, M KK , but far above the masses of all zero modes which would be 
around the weak scale. 

The KK threshold corrections in 5D theories on warped S 1 jZ<i x Z' 2 have been discussed before in [33, 35, 36]. In 
[33, 35], A a for supersymmetric 5D theories on warped S 1 /Z 2 x Z' 2 have been computed in the framework of 4D 
effective supergravity. In this framework, A a could be obtained by computing the tree-level Kahler potential and 
also the one-loop correction to holomorphic gauge kinetic functions which can be determined by the chiral anomaly 
structure of 5D orbifold field theory [37] . The KK threshold corrections for nonsupersymmetric 5D theories on warped 
S 1 /Z2X Z' 2 have been computed in [36] by directly evaluating all KK mode contributions in dimensional regularization 
scheme [38], and it was confirmed that the results in supersymmetric limit agree with those of [33, 35]. In this paper, 
we compute the KK threshold corrections in generic 5D theories on flat S 1 /Z2 x Z 2 with quasi-localized zero modes 
using the method of [36] . 

To calculate one-loop gauge couplings at low energies, we integrate out all massive KK modes and derive the 
one- loop effective action of gauge field zero mode A^: 



S c s = Sb 



(37) 



where Sbarc is given in (34), and the 1-loop contributions from 5D scalar (f>, 5D spinor if and 5D vector Am are given 

by 



iT d 



1 



Tr ln(-C 2 + M 2 (</>)) : 



iT 9 = -Tr*ln(-Z? 2 + Af 2 (*) + ^Jf; 2 )+Tr* ln(^), 

iV A = -^Ti A (- D 2 + M 2 (A^) + jr)-l^ A5 In (-D 2 + M 2 (A 5 j) 
+ Tr^ln(-i} 2 + M 2 (£)). 



(38) 



Here Tr$ implies the functional trace for the 5D field <E>, so contains the summations over the whole KK modes. 
M 2 ($) denotes the mass-square operator whose eigenvalues correspond to the KK mass spectrum of $, and Jj l!/ is 
the 4D Lorentz generator for spin j, normalized as tr( J^ v J^ a ) = C{j){g^ p g ucr — g^ a g up ) where C(j) = (0, 1,2) for 
j = (0, 1/2, 2). Tr^ denotes the functional trace for the zero mode of \E" and Tr^ is the trace over the ghost fields. 

The most convenient way to calculate Tr$ is to replace the KK summation by a contour integration with an 
appropriately chosen pole function P(q). The pole function we will use here has the form 

N'(q) 



P(Q) 



2N(q) 



(39) 



where N(q) has zeroes at q 2 = M 2 (<!>) — m| for M„($) denoting the n-th KK mass and to$ denoting the bulk mass 
of <f>. Note that our pole function is different from the pole function of [36, 38] as the pole positions are shifted by 
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FIG. 1: Contour Ci in the complex q-plane. Bold dots represent poles at q 2 — M^($) — m% where M n is the KK mass 
eigenvalue and m<s> is the bulk mass of 5D field 3>. Note that a 4D state with M n < m$ appears as a pole at the imaginary 
axis. 



CI 



m|, which is mainly for the simplicity of calculation. Using the analysis of Sec. II, we find the following forms of 
N- functions for 5D scalar and fermion fields: 



JV+(«) = 


-i(q 2 + sin(g7ri?) + 


JV-(g) - 


-(/X-l sin(g7ri?) + qcos(qirR)) , 


JV+(«) - 


- (— sin(<77ri?) + q cos(qnR)) , 


N4-{q) = 


- sinfWi?) , 


N*++(q) = 


- sin(g7ri?) , 


N 9 +-(q) = 


— ( — A4_| sm(qirR) + qcos(qnR)) , 


JV*-+(9) - 


-(A4 y sm(qirR) + qcos(qnR)) , 


^*~(«) - 


— sin qnR , 

q 



(40) 

where the scalar and fermion mass parameters /i zz ', and are defined in (2) and (13). The TV-functions for 

5D vector fields AjjJ are also easily found to be 

N A++{q) = - sm(qnR) , 

NA+-(q) = cos(qnR) , 
Na-+(q) = cos(qnR) , 

JVU-(g) = isin( g7 ri?), (41) 

where the boundary conditions of A z ^ are given by 

Af (-y) = zAf (y) , A'f (-y + nR) = z'A**' (y + nR) , 
Af (-y) = -zAf (y) , Af (-y + nR) = -z'Af (y + nR) . 
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FIG. 2: The contour Ci on the upper half-plane can be deformed into C2 without touching any singularity. The poles at the 
imaginary axis do not overlap with the brance cut unless there is a tachyon state. 
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For the pole function defined as above, it is straightforward to find 



Trln(-£ 2 + M 2 ($)+F^Jf ) 



(42) 



rf 4 fc (jp + fc) 2 + g 2 + mj) - \ {p + 2kY{p + 2ky 



d(j) 



(2tt)4 

2C(j)(p 2 g^ -p»p v ) 



(fc 2 + q 2 + ml) ((p + fc) 2 + q 2 + ml) 



d*k 



d 4 p 



Ga(p)A a J-p)(p 2 g 



2 n \x v 



(2^) 4 (fc 2 + q 2 + m 2 ) ((p + fc) 2 + <? 2 + ml) 



(43) 



where d(j) = (1,4,4) and C(j) = (0,1,2) for j = (0,1/2,1), and T a (X) = Tr (T 2 (X)) is the Dynkin index of the 
gauge group representation X. Here the contour C\ is depicted in Fig. 1. To regulate the divergent part of the above 
integral, we split the pole function into two parts: 



P(q) = P(q) + P 00 (g), 



(44) 



where P — > 0(q 2 ) at |g| — > oo. Then is given by 



/4 ? 7T /? 

Pco(9) = r e(Im(9)). (45) 

where e(x) = x/|x| and A is a real constant depending on the Z 2 x Z' 2 parity of the corresponding 5D field: 

A = (-1/2, 0, 0, 1/2) 

for Z 2 x Z 2 parity (Z$,Z'$) = (++, H — , — h, )• With the decomposition (44), all UV divergences appear in the 

contribution from P^ in a manner allowing simple dimensional regularization. 

The 4D momentum integral d 4 k in (43) exhibits a branch cut on the imaginary axis of q for p 2 > 0. For the 
contribution from P, one can change the contour as in Fig. 2 since the contribution from the infinite half-circle 
vanishes. Note that the poles of P on the imaginary axis do not overlap with the branch cut as long as there is no 
tachyon. After integrating by part, we find that the part of Q a from P is given by 



AGa(P) = 



T a ($) 
8tt 2 

r„($) 

8tt 2 



ld(j)-2C(j)) Hi) 



q—noo 



j\x(^d(j)(l-2xf~2C(j)^(q) 



(46) 
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where 

.Ffo) = - m A + Aln g + — g. 

In fact, F{q) turns out to be vanishing in the cases we are now considering. The contribution from P^ includes 

the log divergence from the pole term 1/q. This can be regulated by the standard dimensional regularization of 
4D momentum integral, d 4 p — ► d D p, yielding a 1/(1? — 4) pole. On the other hand, the step- function contribution 
from e(Im(q)) involves a 5D momentum integral which is linearly divergent, but it simply gives a finite result in 
dimensional regularization. Note that linearly divergent correction to l/g% depends highly on the used regularization 
scheme. Namely the coefficient j a of Eq. (36) is regularization scheme dependent, and j a — in the dimensional 
regularization scheme we are currently using. However this does not have any special meaning since the physical 
amplitudes are always expressed in terms of the scheme- independent combination + -^sAttR. Adding the divergent 



9t a 



contribution from P^ to the finite part AQ a from P, we obtain 



Ga 



Ta($) 

8tt 2 



+A 



dx 



1 



1 



dx 



-d(j)(l-2xy + 2C(j) -InN 



-d(j)(l-2xf + 2C(j) 



q=i s /x(l-x)p 2 +m% 



D — A 



(47) 



Using the above result, we find that the KK threshold corrections from a 5D complex scalar 0, 5D Dirac fermion 
\I> and 5D vector Am are given by 



-InA 
P 



T a (<j) ++ ) | In A- 3 J duF(u) In N, 
-3T a (<f>+~) j\uF(u)lnN^+- ^ 
-3T O (0"+) j\uF(u)lnN^- + t 



,++ i 



u 2 p 2 „ 
— + m ++ 



2 2 

— +m\_ 



u 2 p 2 

— + m -+ 



A 



A a (*)+6*ln- = 



jlnA + 3^ duF(u)\nN, 
T a (*++) |-21np + 3^ d«G(u)InJV, 
+ T a (*+-) j+3^ duG(u)lnNy+- ^ 
+ T a (*-+)|+3^ duG(u)lnN^-+ 
+ T a (*~) j-21np + 3^ duG(ti)lnJV*— ^ 



— 1- m_ 



+ M 2 _ + 



u 2 p 2 



M 2 



A a (A M ) + &fln- 
V 



1 

12 



T a (A++W -231nA + 441np + 



/ dwA( 



u)\nN A++ [ —yjp 2 



+ T a (A + M ) f 1 duK(u) In N A+ - 
Ja 



vu, 



T a (A M +) / d«A-(u)lniV A - + (-v^ 



vu 



+ T a (A M ~) |231nA — 21np + y duK{u) In A^- 

where 6f , bf and 6^ are the beta function coefficients for the massless modes from <f>, ^ and Am, respectively, and 

F(u) = uil-u 2 ) 1 ' 2 , 
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G{u) = «(1 - u 2 ) 1 ' 2 - u(l - u 2 )~ 1/2 , 
tf(u) - -9w(l-u 2 ) 1/2 + 24u(l-u 2 )- 1 / 2 . 

For the case that is much smaller than the lowest nonzero KK mass, the above results are simplified to yield 
A Q = g [T a (A++) + T a (A M -)] ln(A7ri?) 
-iT a ($++) l n f Me m++ " R - e— 



6 v ' V 2m 



i To (^++) ln f (^++ + - /i / ++ )e m++7rfl - (m ++ - /i++)(m ++ + ^ ++ ) e -">++^ 



6 vr ' V 2m ++ A 
_I To (0+-) i n A m +" + /"+-)e m+ " 7rii + (m+_ - M+-)e- m+ -^ v 



6 \ 2m + _ 

6 a[<P ' n \ 2to_+ 
-~T a (<f>— )ln' 



6 \ 2m 

2 , /A(e 7W ++ 7rii -e-- M ++ 7rii ) 
T a (* ++ )ln' 



3 v 'V 
-|T a (*+-)ln(e-^+-^) 

-^T a (*-+)ln(e^- + ^) 

-|T a (*--)In(^ — - 1 ) , (48) 



3 v ' V 2M— 

where m zz >, fj, zz > and fi' zz , denote the bulk and brane masses of <j>, and A4 ZZ > is the kink mass of \P. Here <J> ++ is a 5D 
complex scalar field having a zero mode, i.e. a scalar field with /U = // = to, and (f> ++ stands for complex scalar fields 
without zero mode. The 4D beta function coefficients b a are given by 

ba = -yT a (A++) + \T a {A- M -) + ±T a (<P++) + |r a (* ++ ) + |t o (*~), (49) 

which can be easily understood by noting that gives a massless 4D vector, a massless real 4D scalar, and 
^ ±=t a massless 4D chiral fermion. Then comparing the above A a to the expression of 4D Yukawa couplings in (28), 
one easily finds that generically 

A a = O (lnfoiw)) (5°) 

if the 4D Yukawa couplings y ipq are generated by quasi-localization. 

It is straightforward to find an expression of A a for supersymmetric 5D theories using the above result. In super- 
symmetric theories, there can be two type of bulk fields: vector multiplet V containing a 5D vector Am, a Dirac-spinor 
A and a real scalar S, and hypermultiplet H containing a Dirac spinor vff and two complex scalars (f>,<j)'- The mass 
parameters and Z 2 x Z' 2 boundary conditions of component fields are given by 

V zz ' = { A Z M ' , X zz ' (M = 0), £ 52 ' fa = p! = m = 0) } , 

H"' = {(j> zz '((i = V' = m = -M), <f/ Si ' (ji = = m = M), V ZZ '(M)}, (51) 

where z = — z = ±1, z' = —z' = ±1. Here to, fi and \j! denote the bulk and brane masses of scalar field, and M. is 
the kink mass of Dirac fermion field. We then find 



(A Q ) SUSY = [T a (V++) + T (V— )] ln(A7ri?) 

-T a( ^ + )in( A ^ + ^;;- M++ ^)) 



-T (W+-)ln(e- 



A4 + _7T.R\ 
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-T (W- + )ln(e 
-T a (H ) In 



2M — 



(52) 



and the 4D beta function coefficients 



(Msusy = -3r a (V ++ ) + T (V~) + T a (W++) + T a (W~). 



(53) 



In fact, one can obtain (A a )suSY using the 4D effective supergravity method discussed in [33, 35]. We confirmed that 
the result from 4D effective supergravity agrees with the above result which was obtained from a direct calculation 
of KK threshold correction. 

With the above result on A a , one can obtain the low energy gauge couplings at p < Mkk which are determined as 
(36) at one-loop approximation. In most of 5D orbifold field theories, we have 



M K 



K 



M w 



> 



A 



M KK 



(54) 



by many orders of magnitude, where Mw is the weak scale and Mkk is the lightest KK mass. Then the dominant 
part of higher order corrections (beyond one-loop) to low energy couplings at Mw come from the energy scales below 
Mkk, which can be systematically computed within 4D effective theory. To include those higher order corrections, 
one can start with the matching condition at Mkk- 



1 



gl{M KK ) \g, 



1 



a / bare 



l 

8^ 



A a + b a In 



A 



M 



KK 



(55) 



where A a are given by (48), and then subsequently perform two-loop rcnormalization group (RG) analysis over the 
scales between M K k and M W - If there is a massive particle with mass M between M K k and M w , one needs to stop 
at M to integrate out this massive particle, which would yield a new matching condition at M. For a given model, 
one can repeat this procedure to find the gauge couplings at M w , and compare the results with the experimentally 
measured values. 



In some case, there can be another large mass gap between the lightest KK mass Mkk and the next lightest KK 
mass M' KK . For instance, as we have noticed in Sec. II, has two light KK modes with a Dirac mass 2Aie~ M7TR 
when its kink mass MnR ^> 1. Then the lightest KK mass is given by Mkk = 2Me~ MnR , while the next lightest 
KK mass M' KK = 1/R which corresponds to the mass of the first KK mode of gauge fields. If A4ttR is large enough, 
so that 



^MttR 



M'kk 
M KK 2MR 



> 



A 



M 



= AR, 



(56) 



KK 



the next important higher order corrections would come from energy scales between Mkk and M' KK . Those next 
important higher order corrections can be included by performing the two-loop RG analysis starting from M' KK . The 
corresponding matching condition at M' KK is given by 



9l{M' KK ) 



\ J L 

9a J bare 



8^ 2 



A' + b' In 



A 



M'kk 



(57) 



where 



b'a 



b a + Sb a , 
A a - 8b a In 



A 



M 



KK 



(58) 



for A a given by (48). Here b a denote the one-loop beta function coefficients for zero modes, while Sb a denote the 
coefficient for the lightest KK states. 
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V. APPLICATION TO ORBIFOLD GUT 



In the previous section, we have discussed one-loop gauge couplings in generic 5D orbifold field theory, including 
the KK threshold correction A a as 



1 / 1 \ 1 



A a + b a In - 
V 



(a =1,2,3). (59) 



9aip) V 9a /bare 

The bare couplings here consist of several pieces which are not calculable within orbifold field theory: 



2a 7 bare" <?5a 24 ^ 



2 +TT7^ A ^ + ^ + <- ( 6 °) 



So although it is an well-defined relation between the bare parameters and measurable quantities, (59) does not give 
any useful prediction unless additional information on bare couplings is provided. In orbifold GUT which is strongly 
coupled at A, g\ a and j a are universal as a consequence of unified gauge symmetry in bulk, and both n a and n' a are 
of the order of 1/87T 2 as the theory is strongly coupled at A [29, 30]. We then have 

£L-isb +0 (s*) <~ w >- (61) 

With this information on bare couplings, one would be able to predict for instance the value of QCD coupling constant 
at Mz in terms of the measured values of the electroweak coupling constants at Mz and the KK threshold corrections 
A a computed in the previous section. 

Let us discuss in more detail the effect of KK threshold on the predicted value of the QCD coupling constant. To 
this end, it is convenient to consider 

E T~T . (62) 
where a a = g^/Air and n a are the coefficients determined by 

773 = 1, E^ =0 ' & a6a = - ( 63 ) 



It is then straightforward to find 



a 3 (M z ) \a 3 (M z )J 2ir 
where 



-!- (771 Ai + n 2 A 2 + A 3 ) (64) 



J LE 



K a 3 (M z )J \ai(M z ) a 2 (M z ) 
for the experimentally measured electroweak coupling constants a.\{Mz) — 0.0169 and a 2 {Mz) = 0.0338, and 



(65) 



JLE 



? a a (M KK ) S a a ( Mz ) (66) 



can be determined by the RG analysis below Mkk once the zero mode spectrums are known. Note that when 
expanded in powers of A a /87r 2 , (5le is independent of A a at leading order. If zero mode spectrums correspond to the 
standard model (SM), we have 

115 333 
111 = 218' ^-218 (67) 



and 



a 3 {M z )J 



15.3 + 01 -) (68) 
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for Mkk — 10 13 — 10 15 GeV. Here we consider a rather wide range of Mkk to cover the case that the lightest KK 
mass is suppressed by a small localization factor as Mkk = 2M.e~ MlxR . It turns out that the numerical result is 
insensitive to Mkk, e.g. the variation of 1/03 is within 0(1/tt) even when Mkk varies by several orders of magnitude. 
In the above, we have used the approximation scheme to include two- loop RG evolution below Mkk, and then the 
uncertainty of 0(1/ V) is from the dependence of (5le on Mkk as well as from the piece higher order in A a /87r 2 . For 
more interesting case that zero mode spectrums correspond to the minimal supersymmctric standard model (MSSM) , 

Vi = y, V2 = -— (69) 



and 



7.8 + (-) . (70) 



Here we have assumed the superparticle masses Msusy at 0.3 ~ 1 TeV, and then the uncertainty of 0(1/ it) is mainly 
from the variation of superparticle threshold effects at Msusy- 

To see the importance of KK threshold corrections more explicitly let us consider a class of 5D SU (5) orbifold 
GUTs whose effective 4D theory is given by the MSSM. To break SU(5) by orbifolding, Z 2 x Z' 2 is embedded into 
SU(5) as 

Z 2 = diag(+l,+l,+l,+l,+l), 

Z 2 = diag(+1,+1, +1,-1,-1), (71) 
loading to the following orbifold boundary conditions of the 5D vector multiplct V containing the SU (5) gauge fields: 

V = (8, l)( ++) + (1, 3)( ++) + (1, l)( ++) + (3, 2)%-jl + (3, 2)<+-> , (72) 

where the SU (5) adjoint representation is decomposed into the representations of the SM gauge group. Obviously 
then the bulk SU(5) is broken down to SU(3) x SU(2) x U(l) at y = ttR, while it is unbroken at y = 0. 

The model contains matter hypermultiplets ^ P (5), F p (5), T p (10) and T p (10) (p = 1,2,3) with kink masses Mf p , 
M-f^, M-t v and M.t> v , and also the Higgs hypermultiplets H(5) and H'(5) with kink masses M. H and A4h', where 
the numbers in bracket mean the SU(5) representation. We assign the Z 2 x Z' 2 parities of these hypermultiplets as 

Z 2 (F p ) = Z 2 (F' p ) = Z 2 (T p ) = Z 2 (T' p ) = Z 2 (H) = Z 2 (H') = 1 

Z' 2 (F p ) = -Z' 2 (F' p ) = Z' 2 (T p ) = -Z' 2 (T' p ) = Z' 2 (H) = Z' 2 (H<) = -1 , (73) 
and then the orbifold boundary conditions of matter and Higgs hypermultiplets are given by 

\(+-) , n 



F = (3,lY 1 }-> + (1,2)^ 
F' = (3, l)<+ +) + (1,2)^1 



T - 



(3,2)(+ +) + (3,l) ( _ + 2 7i + (l,l)l + - ) , 
T' = (3,2)<; ( r ) + (3,l)L + 2 ;i + (l,l)l ++) , 
^=(3,l) ( _ + 1 7i + (l,2)(+ +) , 

H> = (3, !)(+-> + (1,2)^. (74) 



Then using (70) and (52), we find 



7.8 + -L [ A gaugc + A higgs + A mattor ] +0[-\ (75) 



a 3 (M z ) 2tt L sauBe —J ~ y n 

where 



— A A 



Agauge = In(TriiA) (76) 



15 



corresponds to the KK threshold correction from the 5D vector multiplet, 



A 



9 

14 



In 



/ sinh ttRMh 



V ttRMh 

is the correction from Higgs hypermultiplets, and 



ttRMh + In 



/ sinh ttRMh' 
V ttRMh' 



ttRMh' 



(77) 



^matter 
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E 



In 



In 



sinh ttRM p p 
ttRMf p 

sinh ttRMt p 
itRMtt, 



(sinh irRM f 1 \ 
— '- - ttRMf' 
ttRMf; J 

( sinh ttRMt' \ 

irRM Tp - In — - ttRMt> 

p \ ttRMt' / p 



(78) 



is the correction from matter hypermultiplets. 
For A7ri? w 10 2 , we have 



2 7r ^s au s c ~ 0-8 ■ 



Then the predicted value of l/a 3 would be very close to the experimental value: 

1 



a 3 {M z ) 



= 8.55 ±0.15. 



(79) 



(80) 



cxp 



if Ahiggs + A mat ter is negligible. In other words, the KK threshold correction from Higgs and matter hypermultiplets 



should be negligible, i.e. 



2^ ( Ahiggs + A ma ttor ) ^ O ^-^ , 



(81) 



in order for the 5D orbifold GUT under consideration to be consistent with observation. 



Obviously, for the class of models under consideration, a hypermultiplet with MttR 3> 1 gives a large threshold 
correction A(l/a3) = 0(A4t:R/2t:) which can make the model inconsistent with the observation. For instance, in the 
model of [9], the Higgs zero modes are quasi-localized at y = by having Mh^R = 11.5 and Mh'kR = 6.9, the 1st 
and 3rd generation matters are brane fields at y = irR and y = 0, respectively, and the 2nd generation matters come 
from hypermultiplets with vanishing kink masses. This model then gives l/as(Mz) — 10.9 which is too large to be 
consistent with the experimental value. 



A simple way to avoid a too large Ahi ggs is to assume that both of the Higgs hypermultiplets have MwR <C 
for instance Mh^R w Mh'kR ~ —10. In this case, the Higgs zero modes are localized at y = ttR, and 



1 



A 



higgs 



9 
147 



In (^JMhMh> kR) w -0.45, 



(82) 



which is small enough not to spoil the successful gauge unification. However for matter hypermultiplets, to generate 
the hierarchical 4D Yukawa couplings through dynamical quasi-localization, one needs to localize heavy and light 
generations at different locations. This means that some kink masses should be positive, while some others are 
negative. One then needs a nontrivial cancellation between the corrections from different matter hypermultiplets in 
order for A mat tor ~ 0. In this regard, an interesting possibility is that 



Mf„ — Mf> , Mr. = Mi 



(83) 



which obviously lead to 



2tt 



(84) 



The relation (83) between hypermultiplet masses can be considered as a consequence of global SU(2)h symmetry 
under which (F p ,Fp) and (T p ,Tp) transform as a doublet. This SU(2) H symmetry is broken down to U(1) 3 b+l a t 
y = ttR by the orbifolding boundary conditions (73). It is then straightforward to introduce a dynamics at y = ttR 
which breaks U(1)3b+l spontaneously down to the matter-parity (— l) 3B+i . This is in fact necessary to generate 
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TABLE I: A set of hypermultiplet masses which give realistic fermion masses and CKM mixing while satisfying (81) for 
successful gauge unification. Here n = 0, 1, 2. 



M Tv /M 


M T > /M 

V 


M Fp /Mo 


M F i /M 

V 


(4,2,0) 


(4,2,0) 


(n,n,n) 


(n,n,n) 


(5,1,0) 


(3,3,0) 


(2,0,1) 


(1,2,0) 


(3,3,0) 


(5,1,0) 


(0,2,1) 


(2,0,1) 


(6,6,0) 


(2,4,0) 


(4,0,2) 


(0,4,2) 


(5,1,0) 


(3,3,0) 


(3,1,2) 


(1,3,2) 


(3,3,0) 


(3,3,0) 


(1,3,2) 


(3,1,2) 


(2,4,0) 


(6,0,0) 


(0,4,2) 


(4,0,2) 



nonzero Majorana masses of neutrinos. Then the brane interactions at y = irR are constrained only by the SM gauge 
group, 4D N = 1 supersymmetry and the i?-parity (_i) 3S + L + 2s ( s — spin). 

The KK threshold corrections (77) and (78) to l/a 3 from Higgs and matter hypermultiplets have a correlation 
with the 4D Yukawa couplings given by (28). So the condition (81) for successful gauge unification provides some 
restriction on the possible forms of 4D Yukawa couplings. However still it is not so difficult to construct models to 
produce the correct form of Yukawa couplings through quasi-localization, while satisfying (81). To see this, consider a 
class of models with Mh^R <C —1 and Mh'kR <C — 1, in which the Higgs zero modes are quasi-localized at y = ttR. 
The quark and lepton Yukawa couplings are assumed to arise from the following brane interactions at y = ttR: 

J d 5 x J d 2 9 5(y - ttR)^ [\ U pq HQ p U c q + \°H'Q p D c q + \* q H>L p E c q } , (85) 

where the Higgs doublets H\ and H2 come from H and H', respectively, the lepton doublets L p are from F p , the 
lepton singlets are from T' p , the quark doublets Q p are from T p , the quark singlets U p and D p are from T' and 
F^, respectively. Then according to the discussion of Sec. II, the physical Yukawa couplings of quarks and leptons 
are given by 

y u pq = yJz(MH)Z(M Tl )Z(MT q )X^ , 
y D pq = ^/z(M H ')Z(M Ti )Z(M F ,)\° , 

V E pq - ^1 Z(M H >)Z(M Fp )Z{M Tq ) Ap 9 , (86) 

where 

z( ^ ) = g l W*-r (87) 

If we further assume that all hypermultiplet masses are quantized in an appropriate unit, then the hypermultiplet 
masses tabulated in Table I give the correct quark and lepton masses as well as the correct CKM mixing angles, while 
satisfying (81) for successful gauge unification. Here the unit mass Mo is defined to given by e - M o^ R - = Cabbibo 
angle 0.2. 



VI. CONCLUSION 



In this paper, we have examined the KK threshold corrections to low energy gauge couplings g\ from bulk matter 
fields whose zero modes are dynamically quasi-localized to generate hierarchical 4D Yukawa couplings. We derived the 
explicit form of threshold corrections in generic 5D orbifold field theory on S 1 /Z2 x Z' 2 , and found that the corrections to 
1 /g\ are of the order of \n(y) /8n 2 where y denotes 4D Yukawa couplings generated by quasi-localization. So generically 
quasi-localization significantly affects gauge coupling unification. We then applied the results to 5D orbifold GUT, 
and discussed the conditions for a 5D GUT to generate hierarchical Yukawa couplings without spoiling successful 
gauge coupling unification. Some examples of such 5D GUTs are presented in Table I. 
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